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Abstract 



A quantum interference mechanism of the stripe phase instability in 
quasi one-dimensional (ID) repulsive electron system is proposed. The 
leading spin-charge coupling term in Landau functional is derived mi- 
croscopically. It is shown that away from half filling, periodic lattice 
potential causes cooperative condensation of the spin and charge super- 
lattices, which constitute a ID analogue of the experimentally observed 
quasi 2D structures in lanthanum cuprates and nickelates. Renormal- 
ization group analysis qualitatively supports the Hartree-Fock picture. 
PACS numbers: 64.60.-i, 71.27.+a, 74.72. -h, 75.10.-b 

Stripe phases recently observed in doped antiferromagnets [|J attract attention 
to the problem of multi-mode instabilities in the interacting electron systems |2[- 
HJj. Spin and charge modulations found in the doped layered transition metal oxides 
possess locked mutual phase and doping dependent periods, incommensurate with 
the period of the underlying crystal lattice. Period of the spin modulation is twice 
the period of the charge modulation. An onset of this picture in CuO"2 and MO2 
atomic planes is detected by neutron scattering experiments [[|] which, below some 
temperature, reveal two systems of intense superlattice peaks. Locations of the spin 
and charge peaks in the two-dimensional square Brillouin zone are: (1/2 ± e, 1/2) 
for spin and (±2e, 0) for charge; as well as at tt/2 rotated positions: (1/2, 1/2 ± e) 
and (0, ±2e). Here units of 2n/a are used, and a is the lattice period. For small 
enough doping concentration, x^, a relation holds: e ~ x^, and at greater x d 's a 
shoulder developes. Commensurate antiferromagnetic order at half filling (i.e. at 
Xd = 0) is characterized with the wave vector Qaf — (1/2, 1/2). Existing theoretical 
works on the stripe phase instability either find variationally the ground state of 
the microscopic models |2|||, or use phenomenological Landau- Ginzburg functional 
constructed by symmetry considerations Q. 

In this Letter we propose a quantum interference mechanism of the stripe phase 
instability in (quasi) one-dimensional (ID) electron system, brought about by the pe- 
riodic lattice potential. A microscopic derivation is presented of the main spin-charge 
coupling term in the Landau functional, which was phenomenologically introduced 
in [|J. A model considered below was already investigated in the Hartree-Fock ap- 
proximation || , but proposed here mechanism of the stripe phase ordering was not 
discussed. Neither was it discussed in the previous studies of the ID electron systems 
within one- loop renormalization group approach || ("parquet" formalism), which is 
considered in the last part of this Letter. 

We start from the Hartree-Fock approximation. The Hubbard Hamiltonian with 
the hopping integral t and on-site repulsion U (> 0) may be written in the form: 
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H = tYl cU> + (\n 2 i ~ {St) 2 ) (1) 

Here an identity: nj-j-n^ = ^n 2 — (5!f) 2 was used. Operators n = n-f + nj and are 
fermionic density and spin (z-component) operators respectively, and a is spin index. 
Hamiltonian ([[]) has convenient form for the Hartree-Fock decoupling in the presence 
of two, i.e. spin and charge, order parameters. At half filling, i.e. at one electron per 
site, the system orders antiferromagnetically |||7]]: S* = (m Q /2) cos (Qaf^)- Away 
from half filling, mean-field approximation |5| predicts a spin density wave (SDW) 
instability in the continuous quasi ID system: S z (x) = (m /2) cos (Q-X — 0), with 
the incommensurate wave vector Q_ = 2kp. In the weak coupling limit, U <C t, the 
SDW instability is brought about by the existence of the nested parts of the Fermi 
surface separated by momentum 2kp. The difference Qaf — Q- = 1/2 — 2k F = e is 
expressed via the hole (electron) doping concentration xj > (xj < 0) as: e = Xd/2. 
Effect of the periodic lattice potential away from half filling has been so far neglected, 
both in the mean-field and in the "parquet" approaches 

Consider an incommensurate SDW expression, S z (xi) = (m /2) cos (Q-Xi — (f)), 
on a chain. Here Xi is the "site coordinate" along the chain, and is a phase 
shift. Allowing for the periodicity of the lattice sites, Xi = iV^a; iVj = 0,±1,±2, 
we may rewrite the cosine in the SDW expression in the equivalent form (xj is 
taken in units of a): cos (Q-Xi — 4>) = (l/2)(cos (Q-Xi — (p) + cos (Q+Xi + </>)) = 
cos (xj/2) cos (exi + 0); where Q + = 2ir/a — Q-, i.e. 1/2 + e in our dimensionless 
units. These simple relations demonstrate significant physical fact that an incom- 
mensurate SDW on the lattice can be decomposed into two umklapp related SDWs, 
which could be called "direct" (Q_ = 2kp) and "shadow" (Q + = 2ii/a — 2kp) waves 
P . It is wellknown |10| that formation of SDW with wave vector Q may be regarded 



as Bose-condensation of electron-hole pairs Ck, a c\ + Q a \0) with momentum Q (\0) is 
unperturbed vacuum state of the Fermi-system). In our case this means formation 
of the two electron-hole condensates, with wave vectors Q±, corresponding to "di- 
rect" and "shadow" SDWs. Then, scattering of electrons (holes) by some "extra" 
periodic potential with "matching" wave vector Q + — Q_ = 2e would cause quantum 
interference between the wave functions of the Q± condensates. At some value of 
the relative phase shift the interference may become constructive, causing enhance- 
ment of the gap in the fermionic quasi-particle spectrum of the quasi ID electrons 
and related decrease of their total (free) energy. Here we consider a charge density 
wave (CDW): p(x) = p Q cos (Q + — Q-)x, as an option for the "matching" extra po- 
tential. Depending on the balance between the gain of energy and the "cost" of the 
CDW formation, a condensation of the coupled SDW-CDW phase takes place as a 
first/second order phase transition. This state, consisting of "direct" and "shadow" 
SDWs and "matching" CDW, with locked relative phase shift <f), is right a weak 
coupling analogue of the stripe phase order considered in the literature 0-0- 

In the presence of the SDW and CDW condensates, i.e. densities m(x) and 
p(x) introduced above, the single-particle eigenstates of the Hamiltonian Eq.(|l]) in 
the Hartree-Fock approximation can be determined from the Bogoliubov-de Gennes 
equations derived in pj: 
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=F z2 *-^T + ^"P( X ) M ± _ — ^(a;)M T = £w± (2) 

Here left- and right-movers representation is used for the quasi-particle wave function: 
i/j a (x) = u+(x) exp (ix/A) + au-(x) exp (— ix/4), and coordinate x is expressed in 
units of the chain period, a. The spin density is: m(xi)/2 = (m G /2) cos (Q-Xi — 0) = 
m(xj) cos (xi/2), so that: m(x) = (m G /2) cos (ex + 0). Thus, only slowly varying 
functions u±(x), m(x) and p(x) are involved in Eq.(^). Now, instead of dropping 
CDW term, p(x)u±, from Eq.(^) ||, we shall explicitly allow for it by writing wave 
functions u±(x) in the Bloch wave basis of the periodic CDW potential: 
u±(x) = c^exp (ikx =F iU/(2v) J x p(x')dx') rs c k ±e ikx {J (z) =F J%(z) (e i2ex - e~ i2ex )}; 
where v = 2t is the Fermi velocity of electrons, and z = Up /(2ve). Here J n is the 
Bessel function of the integer order n, and the terms of higher orders than n = 1 are 
neglected, provided that Up /(2ve) < 1. After substitution of u±(x) in the above 
representation into Eq. (0) one finds an algebraic system of linear homogeneous 
equations for the coefficients c^_. Solving corresponding determinant equation one 
finds the single-particle spectrum: 

^ = -y±^(fc+0% 2 + A2; (3) 

= Z7m /?7p \ 
4 J V 2w J 

f{zf = J 2 (z) -2cos(2 ( j ) )Jo(z)J 1 (z) + Jl{z) (4) 

In the "electron doping" case the sign in front of e in Eq.(|3p and of cos 20 in Eq. (|J) 
should be changed. The physical implication of Eqs. @, @) is remarkable. Namely, 
in the presence of the CDW the effective coupling constant U f(z), responsible for the 
SDW condensation, is enhanced with respect to bare coupling U, provided f(z) > 1. 
The latter condition fixes cos(20). The form of Eq.(^j) manifests quantum interference 
between scattering amplitudes of electron in the combined periodic potentials of Q±- 
SDW and (Q+— Q-) CDW. A simple form of solution (^]) is valid in the weak coupling 
limit, U <C t, not too close to half filling, i.e. when Xd 3> A/t. Also, despite the 
gap in the spectrum (||), the system may remain conductive due to sliding of the 
incommensurate density waves along the lattice. 

Free energy of the system (per unit of length), Q, at finite temperature T(= 
follows from Eq.(|3|) and the Hartree-Fock form of the Hamiltonian (J!]): 

n = (U/8){m 2 /2 + pi) - (AT/nv) In [2 cosh (/3£(0/2)] d£ (5) 

Jo 

where E(£) = VF+A 1 , and E b (~ t) is the upper cutoff of the electron energy. Due 
to interference term —2cos(2(p)Ji(Up /2ve)Jo in Eq.([|), the lowest order expansion 
of Q, Eq.(|]), in powers of the CDW amplitude (at small m ) yields linear in p 
term: 5Q ~ —m 2 p U 2 /(tXd)- This term reveals possible microscopic origin of the 
corresponding spin-charge coupling term in the phenomenological Landau- Ginzburg 
functional considered in M. 

Here we merely list the main results following from Eq.(^|). 
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i) Coming from the high temperature limit, A = 0, the stripe phase condenses 
first with cos (20) = —1 or cos (20) = 1 depending on the sign of Xd- Thus, the 
spin density behaves as: (2S z (x)) = — m cos (x/2) sin (ex) in the case Xd > 0, and 
as : (2S z (x)) = m cos (x/2) cos (ex) in the Xd < case. Simultaneously, electron 
charge density is the same in both cases: — p Q cos (2ex) . Hence, the nodes of the spin 
density coincide with the minima (maxima) of the charge density p(x) in the case of 
the hole (electron) doping. This topology is in accord with the stripe phase topology 
considered in the strong coupling limit 0- ||]. 

ii) When doping concentration Xd decreases below x Q ~ yJt/U exp (—2irb/U) the 
system enters strong coupling regime, see Figs. 1,2. This transition is governed 
by dimensionless parameter :U ' p j{2ntxd)- At \J p j [2-ntXd) > 1 higher order Bessel 
functions contribute to the wave function u±(x) and a few-harmonic approximation, 
resulted in Eqs.(|3|), (f|), fails. This indicates a solitonic, rather than SDW-CDW struc- 
ture of the stripe phase in the lowest doping limit, Xd — > 0, in a qualitative accord 
with the previous proposals Stripe phase condensation temperature, T c , mono- 
tonically decreases from the highest value T m = 2(7/7r)t exp (— 27r£/({7/^J) at small 
doping concentrations, \xd\ < x Q , to the lowest value T$dw — 2(j/n)texp (—27it/U) 
at \xd\ ^> x Q . Here 7 = 1.78, and f m = f(z m ) « 1.2 is the maximum value of the 
function f(z) in Eq.([|), reached at z = z m « 0.83. The drop of T c is accompanied 
by a substantial decrease of the SDW and CDW amplitudes at zero temperature, see 
Fig. 1: 

m = 8i/[/exp (— 2irt/U), at x d x Q ; and m = 8t/(Uf m ) exp (— 27rt/(E//^J), at 
x d < x Q . 

Simultaneously, the character of the phase transition changes at x Q from the first order 
(xd < x ) to the second order, Fig. 2. The jumps of the CDW and SDW amplitudes 
at the first order transition temperature are: ml w XdZ m T m \ / 2TrtU / (U f m ) 2 and p Q ~ 
27TXdZ m t/U. The latter estimate for p Q results in the following value of the "small 
parameter": Up /(2ntxd) ~ 0.83 < 1. Hence, our results in the region Xd < x Q , 
based on the neglect of the higher order SDW/CDW harmonics, might be considered 
as qualitative rather than quantitative. 

iii) In the second order phase transition regime, \xd\ S> x G , the order parameters 
close to T c behave as: p a ~ 3rTg DW /(x c itU), and m w 12y/rTsDw/U; in qualitative 
accord with J| (here r = 1— T/Tsdw)- I n the classification of [|f the phase transitions 
described in ii), iii) above belong to the "spin-charge coupling driven" and "spin 
driven" kinds. 

Important issue for the (quasi) ID systems is the influence of fluctuations. We 
study it within a single- loop renormalization group (RG) scheme, so-called "par- 
quet" approximation [H, which we adjust for the case of the two order parameters 
(SDW/CDW) coupled already on the mean- field level. Conventionally, "parquet" - 
RG equations describe behavior of the two-electron scattering vertices 71 (£), 72(^)5 
and 73(0; accounting for back-, forward- , and umklapp scattering of electrons respec- 
tively close to the Fermi "surface" points: ±kp. The RG variable, £, is the logarithm 
of the infrared cutoff of the energy/momentum transfer. It is involved in the (log- 
arithmically) diverging corrections to the vertices, which are initially defined in the 
Born approximation: gi = 7j(£ = 0). Within "parquet" approach only corrections 
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of the highest power in £ are retained in each order of the perturbation expansion 
in each 7, and then summed to an infinite order. An instability of the electron sys- 
tem is manifested by divergences of the (initially finite) vertices 7 »(0 at some finite 
value £ . In the case of the Hubbard Hamiltonian ([!]) at half filling: g$ = U/(4nt), 
i = 1,2, 3, in the dimensionless units. Away from half filling the umklapp condition: 
Pi + P2 = P3 + Pi =t 27r/a, can not be fulfilled when all the quasi-momenta of elec- 
trons (before and after scattering) are close to the Fermi surface. In conventional 
"parquet" theory || an instability does not occur in repulsive ID system without 
umklapp, i.e. when (71,^2 > 0, g 3 = 0. In the (hole) doped case the deficiency of 
momentum transfer: lnja — Akp = 2e = Xd 7^ 0, provides a "natural" momentum 
cutoff, such that at £ > <^ = \n(l/xd) the growth of | 7 i(0l terminates. Hence, in 
the standard approach U, an instability at £ Q would be only possible if £ Q < i.e. 
in the low doping region: Xd < X\ = exp(— £ Q ). In this region commensurate SDW 
is indistinguishable from incommensurate one within a logarithmic accuracy of " par- 
quet" with respect to the infrared momentum cutoff. Using estimate 0: £ c ~ 2nt/U, 
we find: X\ = exp (—27it/U), which is inside our mean-field "strong coupling" region: 
Xd < Xi < x Q . In order to probe the system for a stripe phase instability in the region: 
Xd > Xi, where incommensurability is well resolved, we modify "parquet" treatment 
by adding ficticious "stripe phase" vertices 7«(0 with infinitesimal "starting" values 
li(£,d) at £ d < £o- The vertices 71,2 (0 describe "umklapp" scattering with the wave 
vector 2e = Xd, brought by the CDW component of the stripe phase; while 73(C) is 
due to combined lattice-CDW umklapp: 2n/a — 2e. Thus, "enriched" RG-" parquet" 
equations in the interval £d < £ < 00 become: 

7a(0 = "273(074(0; 74(0 = -4^(73(073(0) 

74(0 = -273(073(0; U0 = -273(074(0 (6) 

where 74(0 = 7i(0 — 272(0; an d same relation is valid between 74 and 
71,2 Q- Diverging solutions of Eqs.(^) for two concentrations Xi > X\ 
are shown in Fig. 3. Analytical solution in the interval £ > £d is: 
73 (0 = J Bcosh(Cln|tanhD(£-£ o )|+0 o ); = ±^2( 7 |_- 5 2 ); T4 (0 = 

(D/2) coth2D(£ - i ); 73(0 = DC/{s/2smk2D{i - Q), where | ~ £ + 0.5(£ o - 
£d) In (2/|7i(£d)|); an d all the constants are determined from the boundary conditions 
for 7j and 7, at £ = and £ ~ £d respectively. 

The divergence of the probe vertices 7i(0 signals 0] in favour of spontaneous 
incommensurate (stripe) ordering in the ground state of the system away from half 
filling. Though, neither mean-field, nor "parquet" approximation gives decisive an- 
swer about the long-range order and/or gap in the quasi-particle spectrum of ID 
system p. 

In conclusion, a quantum interference mechanism of the stripe phase ordering in 
repulsive (quasi) ID electron system is proposed. Preliminary ID renormalization 
group ("parquet") analysis does not contradict the mean-field construction. Sub- 
stantial enhancement of the stripe ordering temperature due to interference between 
charge and spin order amplitudes is found. Analysis of analogous possibility for the 
superconducting order parameter is now in progress. 
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I. FIGURES 



Fig.l. SDW (m ) and CDW (p Q ) stripe phase order parameters, as functions of doping 
concentration Xd at zero temperature, calculated using Eq.([|). Inset: normalized 
stripe phase transition temperature T c as function of doping. 4t/U=3.2. 

Fig.2. Calculated temperature dependences of the stripe phase order parameters, 
m (curves labeled with "m") and p Q (curves not labeled) for the different doping 
concentrations x. Each pair of lines of the same type (e.g. solid, dashed, etc.) show 
m and p Q for each particular value of doping concentration x. 4t/U=3.2 . 

Fig.3. Solid lines: numerical solutions of the Eqs.(|5]) for the vertices 7 3j4 (curves 3 and 
4 respectively) and 73,4 (curves 3' and 4' respectively), for values of ^ = £1,2 < £ D - 
Dashed lines: same as above, but for ^ > (,o- Starting values: gi.2,3 = —Qi = 0.1; 
7ste) = -74(6) = 0.01. 
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